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Abstract 



A specific combination of s-state Lamb sliift AEl{1si/2) — n^AEL{nsi/2) 
is considered. Its value in calculated both in the hydrogen and deuterium 
atoms for n up to 12. The result inludes all correction which can contribute 
1 kHz and particulary: one-loop self energy and vacuum polarization, two- 
loop contribution. Nuclear finite-size corrections for the isotopic difference 
of the combination are also evaluated. 



1 Introduction 



Recently new experimental results on the hydrogen and deuterium Lamb 
shift have been obtained [Q-Pl and some even higher-precision experiments 
are going to be completed |^. The each measurement data include a 
combination of the Lamb shift of several nsi/2 levels (n = 1 12). To obtain 
correlations between theoretical calculations of the Lamb shift of s-states 
with different values of principal quantum number n or their combinations, 
the specific difference of Lamb shifts 

A{n) = AEl{Isi/2) - n^AEL{ns,/2) (1) 

is considered in this work. 

The values A(?t,) have to be used to find self-consistent values of the 
ground state Lamb shift and the Rydberg constant P, H. 

2 General expression and definitions 

The Lamb shift is defined here as a shift from the value of level energy 

E{nl,) = [f{nj) - 1] - ^^^^ [f{nj) - I]' , (2) 



where 



1 , (^«)' 



V 



n - J - 1 /2 + ^(j + 1/2)2 _ ^za)\ 



is the dimensionless Dirac energy with infinite nuclear mass, is the reduced 
mass, Z is the nuclear charge in units of the proton one, and relativistic units 
in which h = c = 1 and a = are used. 

The general expression for eq.(P has the form ^ 



An =^ -olog-rM 1 + +^a'x 



1 



'4( logn - ij{n + 1) + ^(2)) - ^^^^^^) log ^,+AlJ'{n)+G'J{Za) 
lAZm f , , .X , a^iZafm ^ n 1 

-TIT + 1) - - + ^) } + 

(3) 

where = {djdz) logr(z). 
Here 

• logA;o(ns) is the Bethe logarithm, its value can be found in Refs. 

©US; 



• a logarithmic coefficient in the term a(Zaf'\o<^(Zoi)m was obtained 



in Ref . 



• G^iZd) is one-loop self-energy correction in the order of a{Zafm 
and higher; 

• J^^in) is the a;(Za;)^m-contribution of the vacuum polarization; 

• i?62 is the leading logarithmic two- loop correction coefficient; 

• the three-loop term of the order a^{ZaYm in eq.(||) is known and equal 
to zero; 

• recoil corrections in orders (Za)'^vp? jM and iZoif'vp? jM"^ are equal to 
zero; 

• the nuclear finite-size contributions for the difference are small and 
not included. 

The Bethe logarithm is presented in Table 1. The one- and two-loop 
contributions and recoil corrections are discussed below. The nuclear charge 
distribution correction is considered in appendix B. 
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3 One-loop self-energy contribution 



3.1 Extrapolation over Z 

The values of G^iZa) with n = 2-^5 for the hydrogen atom can be derived 
from the numerical data given in Refs. |]13|, [l^ after extrapolation to Z = 1, 
with including for the difference (|I]) only two terms 

G'n'^iZa) = Al,^{n) + {Za)A',,^{n), (4) 
because the logarithmic coefficient in the order a{Zaym has the form |§] 

and does not contribute to the difference A(n). In case of n = 2 a result from 
Ref. [T^ for Aqq{2) has been also used. 



The input data and the results of the extrapolation are given in Table 2. 

The quoted uncertainties arise from statistical errors of the numerical 
integrations of Refs. |T3|, Q and from estimates of the systematic error of 
the fit of eq.(l). 

3.2 Extrapolation over n 

The extrapolated value of G^^{a) can be rewritten in the form 

Gf(«) = ^G„, (5) 

which is more convenient for the extrapolation over n. This transformation 
is needed because the n-dependence of G^^{a) is not well-behaved. That is 
easy to see from the trivial result 

Gf^(a) = 0. 

After extracting the factor {n — 1) explicitly the dependence of G„ over n 
became quite well-behaved, as well as in case of other coeffitients of eq.(|^). 
Our results for n = Q ^ 12 are presented in Table 3. 
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4 One-loop vacuum polarization 



As it was demonstrated in Ref. [|T6], for a calculation of the contribution to 
the difference A{n) it is sufficient to use the approximate wave function 

= ( i ) (6) 

\ 2m / 

where (p is the Schrodinger wave function, and the effective local potential 



VvpM 



a{Za) 



15 



5(r) 



35 



Sir) 



(7) 



The difference coefficient of the vacuum polarization contribution is found 



to be PI 



AVH = ^ (logn - ^{n + 1) + m + ^ + 



(8) 



Results for n = 2 and n = 4 are in agreement with those from Refs. [17 



and ||T^. The approximation of eqs. |^ leads to the same result as in 



work |19| for np,-states. 



5 Two-loop corrections 



The leading two-loop contribution has order a^(Za)^m log^(Za). It origi- 
nates from two- loop self energy of an electron in the Coulomb field H, a 
general expression of which has the form ^2 



5Ef\nl,) = {nl,\L^{Eni^)Gc{Eni^)T.^{E^i^)\nl,) 



dYi ( E) 

+ {nlj I Sa {Eni^ ) I nlj ) + {nl^ | S i ) | nlj ) {nl^ \ \ nlj ) 



E—E„ 



(9) 



where Er(-E') is the r-loop one-particle- irreducible self-energy operator of an 
electron in the Coulomb field and Gc{E) is the reduced Coulomb Greens 
function. 
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In the Yennie gauge |2y, ^ , in which the photon propagator in momen- 
tum space has the form 



the result arises from only the first item of eq.(P) and the ladder part of the 
second one. Contributions of separated diagrams to different energy levels 
are presented in Table 4. 

The contribution was found in Refs. ^ pB , 
form 



and finally it has the 



16 / n — 1 n — 

B,, = -[\ogn-^Pin) + ^P{l)-—- + ^ 



(10) 



The uncertainty of two-loop correction is due to terms beyond log^ Za and 
it is estimated as half the contribution of the leading logarithmic term. 

The same diagrams contribute in the order a'^{Za)^m\og Za also to the 
decay widths of ns- and np-levels (see for details Ref. [0), which can be 
presented as the imaginary part of the self-energy of an electron. This fact 



is useful for checking [|^. An expecit expression of Gc{0,t^', Eni) found in 
Ref. 



SOU has also been used for that. 



6 Recoil corrections 

The evaluation of pure recoil contributions has recently been completed for 
s-levels in the orders {Za)^m^/M || and {Zaym^/M^ |3|]. 



The results for the term of order {Za)^m'^ /M 

*^=M = ^^^(«"2-^). 131]. (11) 

and 

are in disagreement. In our definitions @ both of them do not contribute 
to the difference A(n). Numerical results for recoil correction in the Is- and 



5 



2s-states obtained in Ref. without any expanssion over (Za) confirm the 
scahng factor (see table 5) and they are in fair agreement with eq. ([TI|) . 
It should be also menshioned that results of Yelkhovsky and Pachucki 
[p5|l for higher-/ levels and in agreement one with the other, and with the 
analytic result for p-state of Ref. ||36| and the numerical ones for 2pi/2 and 
2p3/2 of Refs. ||3^ and respectively. 

The correction in order (Za)'^m^/M^ contributes only for nuclear spin 
/ = or 1, but not for / = 1/2 |3[. The result is 



SE{I = 0,1) 



1 (Za) 



-6, 



10- 



(13) 



2M2 n3 

All recoil corrections in this section are equal to zero for the difference (|l|), 
but no proof on these cancelations is not known without direct calculations. 



7 Results 

The differences A(n) and the contributions to it for the lowest ns-levels of 
hydrogen and deuterium are presented in Table 6. 

The uncertainty of the theoretical expression is a rms sum of the uncalcu- 
lated terms of orders o?'{Zaf'm\o<g{Za) and o?[Zaf'm^ which are estimated 
as half the contribution of the leading logarithmic term and the uncertainty 
of the extrapolation for the one-loop self-energy contribution. The uncer- 
tainty of the isotopic difference 

A^^°(n) = A^""(n) - A^^^(n) 

is neglegible and the nuclear finite-size corrections could be important for it. 
They are considered in appendix B and included in the results in Table 6. 
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A Definitions 



The definition of the Lamb shift which is used in this work is due to a shift 
from the energy level given in eq.(§). It is denoted in this Appendix as 
0{nj). 

The expansion of E{nj) has the form 

E = m± [^^y {Zaf E^,. 

An other definition which is also often used include only all Eof, (i. e. 
the exact Dirac energy with the infinite nuclear mass) terms and Ea2 (the 
Schrodinger energy with reduced mass) and a relativistic recoil term £'14: 





{Zaf m 


1 1" 


\ 2n2 J 


2v? M 


j + 1/2 n_ 



E (nj) = m f{nj) + [mji—m] 

Hence the main exluded items in value E^^ are 

2 3 



• (14) 



E 



16 



n^{2j + 1)3 n^{2j + ly rv>{2j + 1) 



1 



and 



24 



2n3 



+ 



j + 1 /2 2n 



Advantages of the first definition of eq.(0) are: extra Eiq terms do not 
contribute to A(n) and to the fine structure. Also extra i?24 terms do not 
contribute to E{ns) (for nuclear spin = 1/2) and to the difference A{n) for 
any nuclear spin. 

Difference between definitions of eq. (0) and eq. (|T1|) in case of the hydrogen 
atom leads to the shifts presented in Table 7. 



B Nuclear charge distribution correction 

Using the nonrelativistic approximation (the same as in case of vacuum po- 
larisation) the nuclear charge distribution correction to the difference can 
easy be found to be (cf. |T6| ) 
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SA{n) = S{ls) X {Zaf 



il){n + 1) — V'(2) — logn 



(n- l)(n + 9) 
ill? 



where £{ls) is the well-known nuclear charge distribution correction to the 
energy of the Is state 

£{ls) = hzafmUr^). 



One can see that the correction to difference is small and can be ne- 
glected both for hydrogen and deuterium atom, but for isotopic difference. 



They are slightly n-dependent (cf. term A 



VP, 



n 



in table 6) and the contri- 



bution to /S.{n) (n = 2 — 8) can be estimated as —0.049(9) kHz for hydro- 
gen, —0.318(33) kHz for deuterium and as —0.269(32) kHz for the isotopic 
difference. The uncertainty takes into acount both: ri-dependence and the 
disagreements in proton [^, |3^, ^ and deuteron pT| , ^ radii measurements. 
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Level 


loP" kn IT?/) 


Ref. 


Is 


2.9841285558 




11 


2s 


2.8117698931 




111 


3s 


2.7676636125 


|11 


4s 


2.7498118405 




11 


5s 


2.7408237279 




101 


6s 


2.7356642069 




10| 


7s 


2.7324291292 




101 


8s 


2.7302672607 


|10| 


9s 


2.7287511660 


m 


10s 


2.7276469387 




10| 


lis 


2.7268177825 


m 


12s 


2.7261793406 




101 



Table 1: Bethe logarithms. 
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Input data for the extrapolation over Z 


Z 




Gl\Za) 


Gl^{Za) 







0.912" 








5 


0.860 (150)^ 








10 


0.707 (11)^ 


0.601 (38)'''^ 


0.492 (75)''''= 


0.412 (113)''''= 


15 


0.619 (1)'' 








20 


0.530'' 


0.436 (5)'''^ 


0.349 (5)'''^ 


0.286 (9)'''^ 


25 


0.445'' 








30 


0.362'' 


0.284 (2)'''^ 


0.220 (2)'''^ 


0.175 (2)''''= 




Results of the extrapolation over 


Z 


1 


0.89 (2) 


0.75 (17) 


0.62 (21) 


0.53 (27) 



Table 2: Extrapolation over Z. " - from Ref. [15]; '' - from Ref. [14]; " - 
fromRefs. [13,14]. 
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n 






2 


0.89(02) 


3.56(09) 


3 


0.75(17) 


3.38(78) 


4 


0.62(21) 


3.31(114) 


5 


0.53(27) 


3.29(167) 


6 


0.46(28) 


3.30(200) 


7 


0.40(31) 


3.30(250) 


8 


0.36(33) 


3.30(300) 


9 


0.33(35) 


3.30(350) 


10 


0.30(36) 


3.30(400) 


11 


0.27(37) 


3.30(450) 


12 


0.25(38) 


3.30(500) 



Table 3: Extrapolation over n. 
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< EiGcEi > 


< > 


E{ns) 







A(n) 


^ (logn - ^{n + 1) + ^(2) - ^) log' 


4 n^-l 1 2 1 d 


E{np) 







T{2p) 





S^^^IXSF I ^^11 3 + 24J' 


r(3s) 


2" In 1 ('91ti 5 J- 387\ /,/( 
55¥^^IZSF l^^^4 + 40 j' 






Table 4: Leading logarithmic contributions of two-loop self-energy in the 
Yennie gauge to different energy levels and decay widths (F). Results are 
done in units of a^{Zafm.'' - Ref. [23]; ^ - Ref. [8]; ^ - Rcf. [24]; - Refs. 
[25,26]; ^ - Ref. [27] (numerically); ^ - Ref. [28] (numerically); ^ - Ref. [29] 
(analytically and numerically);'* - Ref. [30] (analytically). 
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Ref. 


6E{ls) 


6E{2s) 


6E{ls)/86E{2s) 


[ 
[ 


31 
32 
34 




-7.4 
2.8 
-7.1(9) 


-0.93 
0.35 
-0.90(6) 


1 
1 

0.99(13) 



Table 5: Recoil correction in order {Z(xfw?jM in the hydrogen atom in kHz. 
Numerical results of Ref. [34] include also higher order in iZa) corrections. 
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n 




A^^(n) 


A^^(n) 


A^s"^(n) 


A^^"(n) 


A^*°(n) 


2 


39 (1) 


8 


-11(5) 


-187232(5) 


-187225(5) 


7.3 


3 


33 (7) 


8 


-15(7) 


-235079(10) 


-235073(10) 


5.9 


4 


27 (9) 


9 


-17(8) 


-254428(12) 


-254423(12) 


4.7 


5 


22 (11) 


8 


-18(9) 


-264162(15) 


-264158(15) 


4.0 


6 


20 (12) 


8 


-19(9) 


-269747(15) 


-269743(15) 


3.5 


7 


17 (13) 


7 


-19(10) 


-273246(16) 


-273243(16) 


3.2 


8 


16 (14) 


7 


-20(10) 


-275583(17) 


-275580(17) 


3.0 


9 


14 (15) 


7 


-20(10) 


-277221(18) 


-277218(18) 


2.9 


10 


13 (16) 


7 


-20(10) 


-278413(19) 


-278410(19) 


2.7 


11 


12 (16) 


7 


-21(10) 


-279308(19) 


-279305(19) 


2.6 


12 


11 (17) 


7 


-21(10) 


-279996(20) 


-269993(20) 


2.5 



Table 6: Corrections and results for the difference (1) in the hydrogen 
and deuterium atoms and isotopic shift. The one-loop self-energy (A'^(n)), 
the vacuum polarization (A^-^(n)) contribution and the two-loop correction 
A^^(n) are the same in both atoms. 
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dOTTPPtlOTl to 


r?^ • ( AE-^ins) 

IV 1 f 1 J— J I / fcO I 


- AE^Hns) \ fin kHzl 

' * -1 y 1 / J / I ±±j. xvj. J.ZJ 1 


n 


{Zafw?/M 


{Zafrrv^/M^ 


Total 


1 





26 


26 


2 


1 


-13 


-12 


3 


2 


-26 


-24 


4 


2 


-33 


-31 


5 


2 


-36 


-35 


6 


2 


-39 


-37 


7 


2 


-41 


-39 


8 


2 


-42 


-41 


9 


2 


-43 


-42 


10 


2 


-44 


-43 


11 


2 


-45 


-43 


12 


2 


-46 


-44 



Table 7: Recoil contributions to difference of the definitions. 
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